Observation of Strong Quantum Depletion in a Gaseous Bose-Einstein Condensate 
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We studied quantum depletion in a gaseous Bose-Einstein condensate. An optical lattice enhanced 
the atomic interactions and modified the dispersion relation resulting in strong quantum depletion. 
The depleted fraction was directly observed as a diffuse background in the time-of-flight images. 
Bogoliubov theory provides a semi- quantitative description for our observations of depleted fractions 
in excess of 50%. 
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The advent of Bose-Einstein condensates (BEC) in 
1995 extended the study of quantum fluids from liq- 
uid helium to superfluid gases with a 100 million times 
lower density. These gaseous condensates featured rela- 
tively weak interactions and a condensate fraction close 
to 100%, in contrast to liquid helium where the conden- 
sate fraction is only 10% 1]. As a result, the gaseous 
condensates could be quantitatively described by a sin- 
gle macroscopic wave function shared by all atoms which 
is the solution of a non-linear Schrodinger equation. This 
equation provided a mean-field description of collective 
excitations, hydrodynamic expansion, vortices and sound 
propagation |2|. 

The fraction of the many-body wave function which 
cannot be represented by the macroscopic wave function 
is called the quantum depletion. In a homogenous BEC, 
it consists of admixtures of higher momentum states into 
the ground state of the system. The fraction of the quan- 
tum depletion r]o can be calculated through Bogoliubov 



theory: t]q = 1.505^ 



where p is the atomic density 



and a s is the s-wave scattering length 3]. For 23 Na at 
a typical density of 10 14 cm -3 , the quantum depletion is 
0.2%. 

For the last decade, it has been a major goal of the 
field to map out the transition from gaseous condensates 
to liquid helium. Beyond-mean-field effects of a few per- 
cent were identified in the temperature dependence of 
collective excitations in a condensate 0, . The quan- 
tum depletion increases for higher densities - however, at 
densities approaching 10 15 cm -3 the lifetime of the gas 
is dramatically shortened by three-body collisions. At- 
tempts to increase the scattering length near a Feshbach 
resonance were also limited by losses jy, Q • Recently, sev- 
eral studies of strongly interacting molecular condensates 
were performed |8l|9|,ll0(. In lower dimension systems, 
the effect of interactions is enhanced. Strongly correlated 
systems, which are no longer superfluid, were observed in 
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ID systems [n],[^ , and in optical lattices E3,E1- Quan- 
tum depletion in ID was studied in 0, llfi Il7^. where 
condensation and quantum depletion are finite-size ef- 
fects and disappear in the thermodynamic limit. The 
transition between a three-dimensional quantum gas and 
a quantum liquid has been largely unexplored. 

In this Letter, we report the first quantitative study of 
strong quantum depletion in a superfluid gas. Exposing 
atoms to an optical lattice enhances quantum depletion 
in two ways. First, the lattice increases the local atomic 
density, which enhances the interactions (by up to an or- 
der of magnitude in our experiment), ultimately limited 
by inelastic collisional losses. The second effect of the 
lattice is to modify the dispersion relation T(fc), which is 
simply T(k) = h 2 k 2 /2m for free atoms. The effect of a 
weak lattice can be accounted for by an increased effec- 
tive mass. For a deep lattice, when the width of the first 
band becomes comparable to or smaller than the inter- 
action energy, the full dispersion relation is required for 
a quantitative description. 

In addition to enhancing the quantum depletion, an 
optical lattice also enables its direct observation in time 
of flight. For a harmonic trap, the quantum depletion 
cannot be observed during ballistic expansion in the typ- 
ical Thomas- Fermi regime. Because the mean- field en- 
ergy (divided by K) is much greater than the trap fre- 
quency, the cloud remains locally adiabatic during the 
expansion. The condensate at high density transforms 
adiabatically into a condensate at low density with di- 
minishing quantum depletion. Therefore, the true mo- 
mentum distribution of the trapped condensate includ- 
ing quantum depletion and, for the same reason, phonon 
excitations can only be observed by in situ techniques 
such as Bragg spectroscopy In an optical lattice, 

the confinement frequency at each lattice site far exceeds 
the interaction energy, and the time-of-flight images are 
essentially a snapshot of the frozen-in momentum distri- 
bution at the time of the lattice switch-off, thus allowing 
for a direct observation of the quantum depletion. 

The experiment setup is similar to that of our previ- 
ous work 20]: A 23 Na BEC containing up to 5 x 10 5 
atoms in the \F = l,rriF = — 1) state was loaded into a 



2 



crossed optical dipole trap. The number of condensed 
atoms was controlled through three-body decay in a com- 
pressed trap, after which the trap was relaxed to allow 
further evaporation and re-thermalization. A vertical 
magnetic field gradient was applied to compensate for 
gravity and avoid sagging in the weak trap. The final 
trap frequencies were oo x ^y iZ = 2tt x 60,60,85 Hz. The 
mean Thomas- Fermi radius was ~ 12 /am for 1.7 x 10 5 
atoms. 

The lattice beams were derived from the same single- 
mode infra-red laser at 1064 nm used for the crossed opti- 
cal dipole trap. All five beams were frequency-shifted by 
at least 20 MHz with respect to each other via acousto- 
optical modulator to eliminate cross interference between 
different beams. The three lattice beams had a 1/e 2 - 
waist of ~ 90 fim at the condensate, and were retro- 
reflected to form standing waves. The two horizontal 
beams were orthogonal to each other, while the third 
beam was tilted ~ 20° with respect to the vertical axis 
due to limited optical access. One and two dimensional 
lattices were formed using either one or both of the hor- 
izontal beams. The trap parameters were chosen such 
that during the ramping of the optical lattice potential, 
the overall Thomas-Fermi radii remained approximately 
constant in order to minimize intra-band excitations. All 
the measurements were performed at a peak lattice site 
occupancy number ~ 7, as determined by a tradeoff be- 
tween small three-body losses and good signal-to-noise 
ratio. 

The optical lattice was linearly ramped up to a peak 
potential of 22 ± 2 Er in time r ramp , and then lin- 
early ramped back down at the same speed. This ramp 
sequence was interrupted at various times by a sudden 
switch-off of all lattice and trapping potentials (< 1/is). 
Absorption images were taken after 10 ms time of flight, 
reflecting the momentum distribution of the system at 
the instant of release (Fig. [I). Based on the number of 
atoms remaining in the condensate after the full ramp- 
ing sequence (> 80%), we concluded that r ramp > 1 ms 
satisfies the intra-band adiabaticity condition. In the fol- 
lowing discussion, all measurements were performed for 

Tramp 50 HIS. 

The loss and revival of the interference contrast, as il- 
lustrated in Fig.QI has been associated with the quantum 
phase transition from a superfluid to a Mott-insulator 
PJ,|l4|. The presence of sharp interference peaks in- 
dicates coherence and superfluid behavior, whereas the 
presence of a single broad peak indicates the insulating 
phase. However, as we show in this Letter, even before 
the lattice depth exceeds the critical value for the phase 
transition, a diffuse background gradually emerges as a 
result of quantum depletion. The interference peaks rep- 
resent the population in the zero quasi-momentum state, 
and the diffuse background represents the population in 
the rest of the Brillouin zone. We only account for the 
lowest energy band as the population in higher bands 
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FIG. 1: Interference patterns in time of flight: The ramp- 
ing sequence was interrupted as the lattice is ramped up 
(11, 18, 22 E R ) and back down (11, Er). The time of flight 
was 10 ms; the field of view is 861 /mix 1075 /mi. 



remains negligible for our parameters. 

In the time-of-flight images, we masked off the sharp 
interference peaks and performed a two-dimensional 
gaussian fit on the diffuse background peak. After the 
lattice was fully ramped down, most of the atoms re- 
mained in the condensed fraction while a small fraction 
(up to 20%) were heated and distributed across the first 
Brillouin zone likely due to the technical noise and im- 
perfect adiabaticity of the ramp. A linear interpolation 
was used to subtract this small heating contribution (up 
to 10% at the maximum lattice depth) and obtain the 
quantum depletion iV q d/iV, where iV q d is the number of 
atoms in the diffuse background peak (quantum deple- 
tion) and N the total number. 

We performed this measurement for one, two and three 
dimensional optical lattices, and the main results are 
shown in Fig. [3 The quantum depletion became signifi- 
cant for lattice depths of > 10 En for a three dimensional 
lattice (Er = h 2 kf att /2m, where fci at t = 27r/Ai a tt is the 
lattice wave number). Note that the Mott-insulator tran- 
sition starts to occur only at lattice depths > 16 Er (see 
below). TVqd/iV remained small for one and two dimen- 
sional lattices. 

A theoretical description of quantum depletion can be 
derived from the Bogoliubov theory which is the stan- 
dard theory to describe the ground state properties of a 
weakly interacting system. The population in the (quasi- 



) momentum state k is given by 



•opulation 



T{k) + n U - ^2T(k)n U + T 2 (k) 
2^2T(k)n U + T 2 (k) 



(1) 



where T(k) is the kinetic energy, no is the occupancy 
number [per cubic lattice cell of (Ai at t/2) 3 ] in the zero 
momentum state, and U is the on-site interaction energy 
HE 13 EE • Incidentally, Vk is one of the coefficients 
in the Bogoliubov transformation. The total occupancy 
number n is given by the sum of no and the quantum 
depletion: n = no + v k- For a gi yen density n, the 

quantum depletion can be obtained by using Eq. JQ and 
the appropriate dispersion relation T(k). 
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FIG. 2: Quantum depletion of a 23 Na BEC confined in a one, 
two and three dimensional optical lattice: the data points 
with statistical error bars are measured N^d/N; the three 
thick curves are the theoretical calculation of N^d/N using 
Bogoliubov theory and local density approximation. For com- 
parison, also shown are (thin curves): (i) the (smoothed out) 
Mott-insulator fraction Nmi/N based on a mean-field theory; 
(ii) the calculated quantum depletion for a homogeneous sys- 
tem of per-site occupancy number n — 1 and (iii) n — 7. 



A band structure calculation was performed to obtain 
the on-site interaction U (and also the tunnelling rate J) 
as functions of the lattice depths, shown in Fig. [3] In cal- 
culating U, we use a Wannier density function along the 
dimensions with a lattice, and a uniform density in the 
ones without. J is independent of the lattice wavelength 
or atomic mass for a given lattice depth (all energies mea- 
sured in Er). 

The quantum depletion for a lattice of uniform occu- 
pation is obtained by integrating over the first Brillouin 
zone: n = no + J v\ dk. For a sufficiently deep lattice 
(> 5 Er), the dispersion relation is given by: 



T(q) = 4 J sin2 fe + *Er £ ll 



(2) 



where dimensions 1 through d are assumed to have a 
lattice present and q = kAi at t/47r. 

For an inhomogeneous system such as a harmonically 
confined condensate, we apply the result from the homo- 
geneous system to shells of different occupancy numbers 
using the local density approximation (as the dependence 
of the quantum depletion on the occupancy number is 
slowly- varying) . The calculated quantum depleted frac- 
tions are plotted in Fig. [3 The semi-quantitative agree- 
ment between the observed quantum depletion and the 
simple Bogoliubov theory, even for values around 50 %, 
is the main result of this Letter. The remaining discrep- 
ancies may be due to unaccounted heating, a systematic 
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FIG. 3: On-site interaction U and tunnelling rate J for a 
23 Na BEC in an optical lattice at 1064 nm: UdD is for a 
d- dimensional lattice. In a three dimensional lattice, the su- 
perfluid to Mott-insulator transiti on for oc cupancy number n 
occurs when J n — 6(2n + 1 + 2yn(n + 1)) J equals U3B [see 
Eq. (0]. The horizontal locations of the crossing points where 
Jn — U3B are the critical lattice depths. 



overestimate of the background, and the inadequacies of 
Bogoliubov theory for large quantum depletion. 

In free space, the dispersion curve is a continuous 
parabola. Both the number of populated states and the 
population in each state increases with the atomic den- 
sity p, and the quantum depletion 770 is proportional to 
p 1 / 2 . This still holds for shallow lattices (< 5 Er) when 
the quantum depletion (77) does not saturate the lowest 

/ * m *\3/2 

band: r] = r]o I J (cmf and m are the free space 

mean- field energy and atomic mass; ej^ F and m* are the 
lattice-enhanced mean-field energy and effective mass) 

E3- 

The situation changes for deeper lattices as the inter- 
action energy U becomes comparable to the width of the 
first energy band (approximately 4J). In this regime, the 
quantum depletion starts to saturate the lowest band, 
but the higher bands remain virtually empty due to the 
large band gap. While the population iV q d in the low- 
est band continues to increase with the atomic density 
p = n/(Ai a tt/2) 3 , the quantum depleted fraction N^/N 
actually decreases with p [see Fig.EJwhere the calculated 
quantum depletion is bigger for n = 1 than for n = 7 at 
large lattice depths (> 9 En)]. 

The fact that the observed quantum depletion for one 
or two dimensional lattices remained small provides fur- 
ther evidence for our interpretation of the diffuse back- 
ground as quantum depletion. In the dimension with 
a lattice present, the band width is proportional to the 
tunneling rate J which decreases exponentially with the 
lattice depth. The interaction energy U increases much 
slower with the lattice depth. Therefore the flattening of 
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the dispersion relation contributes more significantly to 
the increased quantum depletion. Since this flattening 
does not occur in the dimension without a lattice, the 
quantum depletion for one or two dimensional lattices 
is expected to increase much slower compared to three 
dimensional lattice, consistent with our observation. 

For a three dimensional lattice, quantum depletion and 
the superfluid to Mott-insulator transitions are two con- 
sequences of admixing higher momenta into the many- 
body wave function. For increasing interactions, the 
ground state wave function is first a depleted superfluid, 
then developes strong correlations and suppressed den- 
sity fluctuations and finally turns into the Mott insulator. 
In our measurements, the insulating regions appear fully 
in the diffuse background together with the quantum de- 
pletion of the superfluid regions. A lower bound for our 
measured N q ^/N is thereofre provided by calculating the 
fraction of atoms in the Mott-insulator phase. According 
to the mean-field theory in a homogenous bosonic lattice 
system [110 Sill, the critical value (U/J) c at which 
the phase transition occurs for occupancy number n is 
given by: 

(U/J) c = z[2n + 1 + 2Vn(n + l))] (3) 

where z is the number of nearest neighbors (z = 2d for 
a d-dimensional cubic lattice). For an inhomogeneous 
system such as a trapped condensate, shells of different 
occupancy numbers enter the Mott-insulator phase from 
outside as the lattice potential is increased and U/J ex- 
ceeds the critical values. 

In our experiment, the peak occupancy number is ~ 7. 
From Fig. 03 the critical lattice depths in a three dimen- 
sional cubic lattice for n = 1,2, ...,7 are 16.3, 18.5, 20.0, 
21.2, 22.1, 22.9, 23.6 E R respectively. The integrated 
Mott-insulator fraction Nmi/N as a function of lattice 
depths is plotted in Fig. [2j Instead of a step function 
with jumps at each critical lattice depth, we use a smooth 
spline curve for Nmi/N, which is more realistic given 
the fluctuations associated with the atom numbers and 
lattice depths. The measured N q ^/N was significantly 
greater than Nmi/N. 

In the case of one and two dimensional lattices, a Mott- 
insulator transition would only occur for lattice depths 
much larger than those in our experiment. Note that 
Eq. (J3J) is not directly applicable as the dimensions with- 
out lattice beams in our system are only loosely confined 
and cannot be considered frozen^. In addition, n in 
Eq. (JSJ) is the number of atoms per lattice site and far 
exceeds 7 for the lower-dimensional lattices. 

In conclusion, we conducted a quantitative study of 
quantum depletion in a gaseous BEC through the ap- 
plication of an optical lattice, and found reasonable 
agreement with a model based on Bogoliubov theory in 
the predominantly superfluid regime. A complementary 
study was recently reported by Gerbier et al. . of 
the non- vanishing visibility of the interference peaks in a 



Mott insulator as a result of the admixture of particle- 
hole states. The two works together give a complete de- 
scription of the ground state in both the superfluid and 
insulating phases. More elaborate theoretical treatments 
for the intermediate case have been presented in Refs 
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